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Abstract 

Studies  of  plasma  jet  interaction  with  cross  flows  has  been  made  for  development  of  the 
concept  of  enhancing  of  combustion  and  ignition  at  application  of  plasma  jets.  A  computer 
simulation  of  a  cross  flow  nitrogen  plasma  nitrogen  jet  from  the  divergent  plasma  generator  in 
nitrogen  has  been  carried  out.  This  model  involves  a  realistic  equation  of  state,  plasma  radiation 
associated  with  plasma  formation  inside  the  plasma  jet  as  well  as  coupling  between  the  internal 
nozzle  with  external  flow  fields.  Our  model  facilitates  the  application  of  the  implicit  free- 
Lagrange  method  for  carrying  out  of  computations,  which  have  been  carried  out  for  jet  injection 
angles  135°  into  the  cross  flow.  These  show  the  complex  and  non-homogeneous  structure  of  the 
flow.  They  represent  the  necessary  step  in  understanding  of  applicability  of  the  divergent  plasma 
jets  for  propulsion  enhancement. 

Introduction 

The  concept  of  plasma  enhancing  combustion  and  ignition  by  applying  weak  ionization  is 
being  investigated  at  a  number  of  laboratories,  recently1. 

Calculations  simulating  gas  discharge  devices  relate  to  composition  of  thermal  plasmas  in 
air,  nitrogen  and  hydrogen  as  a  function  of  gas  temperature2  3  show  that  for  typical  plasma  jet 
temperatures  in  the  range  5000-10,000°  K  the  conditions  are  already  sufficient  for  ignition  of 
flammable  mixtures  and  the  role  of  combustion  enhancing  by  plasma  chemical  mechanisms  or 
by  plasma  dynamic  processes  becomes  of  interest.  Our  analyses  show  that  plasma  jet  - 
supersonic  flow  interaction  processes  can  easily  stimulate  the  fuel  activation. 

Solving  this  problem  requires  considering  interaction  of  the  plasma  jet  with  a  cross  flow. 

In  a  two-dimensional  approach  we  model  a  row  of  plasma  generators  with  the  divergent  nozzle 
interacting  with  the  flow  at  some  angle  (0-180°)  to  the  flow,  see  Fig.l4. 
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In  this  our  simulation  as  in4’5,  the  gasdynamic  equations  are  coupled  with  the  2-D  Maxwell 
equations,  and  the  plasma  dynamics  equations,  accounting  for  ohmic  heating  and  radiative 
plasma  cooling.  A  Lagrange  method  based  on  a  completely  conservative  implicit  difference 
scheme  with  an  adaptive  triangular  mesh6  is  applied.  This  was  effective  in  modeling  the  coupling 
of  the  internal  and  external  flows  inside  and  outside  the  counterflow  jet  and  its  nozzle5.  This 
approach  is  applied  in  this  investigation  of  interaction  of  divergent  plasma  jets  with  cross  flows. 

In  the  present  work  we  consider  the  construction  of  the  plasma  generator  with  the 
divergent  inter  electrode  channel  where  the  minimum  nozzle  diameter  is  3  mm.  The  half  angle  is 
6°.  The  nozzle  length  along  the  axis  is  30  mm.  The  cathode  diameter  is  4.5  mm.  The  half  angle 
of  the  conical  cathode  is  75°.  The  distance  between  the  cathode  to  the  nozzle  is  1 .7  mm.  The 
transmitting  area  in  channels  is  equal  to  the  transmitting  area  in  the  nozzle  critical  cross  section. 
In  the  initial  part  of  the  nozzle  takes  place  the  nitrogen  inflow  to  the  plasma  generator.  The  range 
of  the  plasma  generator  power  is  10-15  kW,  the  gas  consumption  range  is  0.2  g/min-0.2  g/s.  The 
energy  input  to  the  gas  takes  place  mainly  at  a  distance  of  5  mm  before  the  cathode. 

The  divergent  plasma  jet  was  chosen  for  its  prospective  properties  shown  in  ’ 

Plasma  generators  injecting  plasma  iet  to  a  crossflow 

The  model  is  based  on  the  following  system  of  radiating  plasma  electrodynamics  equations 
verified  in  2,3 

dp/dt  +  div(pv)  =  0,  ( 1 ) 

p  dv/d t  =  -  grad  p, 

p  d£/d t  =  -  p  div  v  +  j2/a  -  q, 

rot  E  =  -  dB/dt, 

rot  B  =  paj, 

div  B  =  0, 

j  =  oE, 

where  p,  p,  v,  e  are  the  plasma  (gas)  density,  pressure,  velocity,  and  internal  energy;  a  is  the 
electric  conductance  of  the  gas  and  plasma,  q  is  the  specific  power  of  plasma  thermal  radiation, 
E  is  the  electric  field  intensity,  B  is  the  magnetic  field  induction,  j  is  the  electric  current  density, 
and  po=47t  l0“7  H/m.  The  unsteady  cylindrically  symmetrical  approximation  with  respect  to  the 
gas  and  plasma  rotation  was  applied:  v  =  (vr,  v9,  vz),  B  =  (0,  B,p,  0),  j  =  (jr,  0,  jz),  E  =  (Er,  0,  Ez), 


3/3<p  =  0;  the  subscripts  r,  <p,  z  correspond  to  the  radial  r,  azimuth  <p,  and  axial  coordinate  z.  The 
plasma  equation  of  state  can  be  expressed  as 


p  =  a  pkBT/mm,i,  e = pApi  Y~  1 )],  (2) 

where  y  is  the  adiabatic  exponent,  kB  is  the  Boltzmann  constant,  T  is  the  plasma  temperature, 
rrimoi  is  the  initial  gas  molecular  mass,  and  factor  a  takes  into  account  the  free  species  number 
changes  resulting  from  chemical  reactions  (primarily  dissociation  and  ionization).  Functions  y  = 
y  (p,T),  a  =  a(p,T),  a  =  a  (p,T)  and  q  =  q(p,T)  for  nitrogen  plasmas  were  computed  preliminarily 
and  approximated  on  the  basis  of  data  in  Refs. 2  3  Calculations  first  assumed  a  =  1 . 1  and  y-1.4, 
as  in  4  5  but  changed  in  subsequent  in  later  parametric  variations  of  the  temperature  T. 

The  computation  region  in  included  both  the  regions  inside  the  plasma  generator  and 
outside  the  model.  No-slip  boundary  conditions  on  the  walls  were  assumed  for  plasma  and  gas 
and  d/dr  =  0  was  assumed  on  the  axis  of  symmetry.  The  gas  input  azimuth  and  radial  velocities 
Vqo,  Vro  were  determined  via  gas  density  po  from  the  given  gas  flow  rate  m'  and  areas  of  the 
inflow  tangent  orifices  Fz  and  of  the  circular  input  slot  Fs:  v^o,  =  m'/(po  FT)t  v^=  m'Apo  Fs).  The 
input  gas  temperature  T0  =  300  K.  Conditions  p  =  pf,  vz  =  -  vzf  ,  T  =  Tf  on  the  right-hand 
boundary  defined  the  supersonic  ambient  flow.  Conditions  on  the  left-hand  boundary  and  on  the 
peripheral  boundary  provided  a  free  exit  of  the  fluid. 

The  entire  left-hand  wall  inside  the  plasma  generator  had  electric  potential  <p\ ;  in  the  slot 
where  the  gas  is  injected,  the  condition  on  the  magnetic  field  defined  via  the  total  discharge 
current  I  was  applied:  B0  =  |iol/(2ra'o),  here  r<>  is  the  corresponding  radial  coordinate.  All  the 
remaining  walls  had  electric  potential  .  Details  of  computations  in  case  of  internal  problem  are 
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presented  in  .  To  expedite  the  computations,  a  small  background  conductance  o  <  10  cWx  was 
assumed  to  be  in  the  cold  gas;  here  cw  is  the  highest  value  of  plasma  conductance  in  the 
electric  current  channel.  The  computed  distribution  of  electric  current  in  the  rest  of  the  channel 
proved  to  be  practically  independent  of  the  fluctuations  of  the  plasma  conditions  (see  below). 
This  is  a  result  of  the  nitrogen  weak  plasma  conductance  dependence  on  plasma  pressure  and 
temperature  in  the  parametric  region  of  interest.  This  allowed  us  to  “freeze”  this  distribution  and 
to  avoid  its  expensive  recalculation  during  the  computations.  So  we  transformed  the 
electrodynamic  problem  to  the  problem  of  the  gas  heating  in  the  channel. 

Preliminary  calculations  corresponded  to  the  conditions  of  7  8  when  the  plasma  injection 
was  realized  to  dead  air  at  atmospheric  pressure  at  temperature  300  K.  Initial  density  is 


determined  from  the  energy  of  state  of  the  ideal  gas,  the  velocity  of  inflow  is  determined  by  the 
data  of  the  consumption. 

Special  calculations  were  made  on  the  divergent  plasma  generator  injection  to  a  gas  at 
(0=180°).  They  are  necessary4  to  determine  the  distribution  of  parameters  over  the  jet  cross 
section  at  the  outlet  for  the  solution  of  the  cross  flow  jet  interaction  at  angles  (0<6<180°). 

In  Fig  2-6  one  can  see  the  distribution  of  main  parameters  of  the  divergent  plasma  entering 
the  gas  at  (0=180°)  the  gas  consumption  2  g/s.  pressure  of  the  oncoming  flow  P=0.2  atm.  at  the 
plasma  generator  power  15  kW,  typical  time  in  calculations  10'4  s.  It  is  2-D  solution  of  the 
plasma-gas  flow  problem,  including  the  internal  and  external  problem  as  a  whole.  In  this  case  we 
use  the  information  on  thermal  conductivity  and  viscosity  characteristics,  since  the  temperatures 
in  the  plasma  generator  are  high.  As  one  can  see  the  distribution  of  temperature,  pressure  and 
density  is  non  homogeneous  over  the  nozzle’s  exit.. 

About  10  calculated  points  in  the  outlet  channel  cross  section  of  the  plasma  generator  were 
used  for  coupling  plasma  generator  internal  nozzle  flow  with  the  external  flow  as  in  our 
counterflow  jet  study  5. 

In  solution  of  the  jet-cross  flow  interaction  at  angles  (O<0<18O0)  the  jet  outside  the  plasma 
generator  is  considered  to  be  flat,  it  models  a  row  of  nearly  placed  plasma  generators.  The 
parameters  of  the  jet  at  the  outlet  to  the  cross  flow  are  taken  from  the  calculations  for  (0=1 80°)4 
Such  formulation  allows  to  clarify  characteristics  of  the  jet  interaction  with  the  cross  flow, 
spatial  distribution  of  gas  parameters  distribution  and  peculiarities  of  the  jet  influence  on 
oncoming  flow  (presence  of  shock  waves,  dependence  on  outflow  angles). 

In  this  case  we  do  not  consider  the  influence  of  the  viscosity  and  thermal  conductivity 
outside  the  convergent  nozzle  in  the  cross  flow  since  their  influence  at  typical  gasdynamic  times 
and  external  temperatures  is  small. 

Questions  of  turbulence  are  not  considered  on  this  study  of  investigations.  In  this  studies 
we  did  not  consider  dissociated  and  charged  plasma  properties  since  this  question  has  been 
analyzed  in  5  it  was  shown  there  that  nitrogen  atom  concentrations  are  of  about  1018  cm'3  in 
conditions  of  interest  for  the  combustion. 

Discussion  of  obtained  results 

Results  for  three  cases  have  been  obtained. 

Obtained  results  of  the  divergent  plasma  get  parameters  at  injection  to  dead  air 
qualitatively  agree  with  experimental  data7,8  in  temperatures  distribution  (0.3  -  ~45  kK)  over  the 


cross  section  at  the  same  distance  from  the  electrode  and  the  same  plasma  generator  power 
W~50  kW.  Concentrations  of  electrons  were  comparable  with  those  of 8. 

At  Fig.2-6  we  present  characteristics  of  the  profiled  plasma  jet  at  the  nozzle  exit 
outflowing  at  180°into  nitrogen.  Results  of  calculations  show  that  at  time  moments  (4  -  1 0)- 1 0  4 
quasi-  stationary  situation  is  realized. 

At  Fig.  7-11  we  present  results  for  a  divergent  plasma  jet  coming  out  from  the  nozzle  at 
(0=135°).  The  jet  screens  a  flow  and  a  shock  wave  is  formed  in  oncoming  flow  The  oncoming 
flow  turns  at  the  shock  and  flows  along  the  jet.  A  region  screened  by  the  jet  is  formed  behind  the 
jet.  Temperature  picture  show  appearance  of  a  cluster  of  highly  heated  gas. 

In  general  temperatures  and  flow  velocities  in  the  place  of  jet-flow  are  consistently  higher 
than  in  the  case  of  the  get  from  the  plasma  generator  of 5  (a  plasma  channel  with  the  barrel-type 
cavity). 

Summary  and  Conclusions 

We  presented  the  extemal/intemal  theoretical-computation  analysis  of  the  divergent 
plasma  jet  -  crossflow  interaction  for  conditions  of  already  known  and  of  planning  experiments. 
This  preliminary  analysis  is  necessary  gasdynamic  part  of  plasma-  jet  applicability  for  ignition- 
combustion  investigation  problems. 

Equations  for  plasma  inside  and  outside  a  channel  including  transport,  chemical  and 
ionization  processes  were  applied.  Derived  equations  allowed  application  of  the  implicit  free- 
Lagrange  method  to  carry  out  computations. 

Our  computations  show  the  complex  and  non-homogeneous  structure  of  the  flow  for  the 
structure  of  a  flow  outside  plasma  generator,  which  is  different  for  different  types  of  plasma 
generator  geometries.  Obtained  results  show  that  the  divergent  plasma  generator  could  be  more 
prospectus  for  fuel  ignition  and  combustionin  in  the  supersonic  flow  than  that  considered  in  4'5. 
Our  results  are  in  qualitative  agreement  with  results  of  known  experiments  for  divergent  plasma 
jets  entering  the  dead  gas.  They  facilitate  interpretation  of  new  experiments,  particularly  the 
relative  roles  of  plasma  and  gasdynamic  processes. 
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Fig.l  Plasma  jet  from  a  slot  created  by  a 
row  of  plasma  generators 


T>  30.37 
1 .5E+02 
1 .2E+02 
9.9E-KM 
•.1E+01 

C.6E-KH 

5.4E+01 

4.4E+01 

3.5E-KJ1 

2.9E+01 

2.4E+01 

1. 96+01 

1.6E+01 

1 .36+01 

1 .06+01  ! 

•  4E+00 

e.iE+oo 

s.oe+oo 

4.3E+00 

3.7E+00 

3.0E+00 


z 


Fig.3  Temperature  distribution  at  the 
divergent  plasma  jet  entering  the  incident 
gas,  consumption  of  a  gas  0.2  g/s,  inflowing 
gas  pressure  0.2  atm,  plasma  generator 
power  15  kW,  time  30- 1  O'4  s 


Fig.2  Pressure  distribution  at  the  divergent 
plasma  jet  entering  the  incident  gas, 
consumption  of  a  gas-  0.2  g/s,  inflowing 
gas  pressure-  0.2  atm,  plasma  generator 
power -15  kW,  time -301 0"4  s 
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Fig.4  Mach  number  distribution  at  the 
divergent  plasma  jet  entering  incident  gas  , 
consumption  of  a  gas  -0.2  g/s,  inflowing  gas 
pressure  -0.2  atm,  plasma  generator  power  -15 
kW,  time -301  O'4  s 


Fig.9  Jet  interaction  with  the  cross  flow  at  (6=135°).  Module  of  the  velocity  distribution  at 
the  divergent  plasma  jet  entering  the  incident  gas,  consumption  of  a  gas  is  0.2  g/s,  inflowing 
gas  pressure  is  0.2  atm,  plasma  generator  power  is  15  kW,  time  is  30- 10  4  s 


